Dynamical Critical Properties of the Random Transverse-Field Ising Spin Chain 
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We study the dynamical properties of the random transverse-field Ising chain at criticality using a 
mapping to free fermions, with which we can obtain numerically exact results for system sizes, L, 
as large as 256. The probability distribution of the local imaginary time correlation function 5(r) is 
investigated and found to be simply a function of a = — log 5(r) / log r. This scaling behavior implies 
that the typical correlation function decays algebraically, S'typ(r) ~ r~°'yp, where the exponent Otyp 
is determined from P{q), the distribution of a. The precise value for atyp depends on exactly how 
the "typical" correlation function is defined. The form of P[a) for small a gives a contribution 
to the average correlation function, S'av(T), namely S'av('r) ~ (logr)~^^™, where Xm is the bulk 
magnetization exponent, which was obtained recently in Europhys. Lett. 39, 135 (1997). These 
results represent a type of "multiscaling" different from the well-known "multifractal" behavior. 
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I. INTRODUCTION 

Quantum phase transitions show a number of remark- 
able features and have attracted considerable interest in 
recent years. They occur at zero temperature and are 
driven by quantum, rather than thermal, fluctuations. 
Hence, they are induced by varying a parameter other 
than the temperature, such as an applied transverse mag- 
netic field. In particular, systems with quenched dis- 
order show surprising properties near a quantum criti- 
cal point. For instance, it was found that in quantum 
Ising spin glaasesEro and random transverse-field Ising 
ferromagnetsj3a as well as in "Bose glass" systems (which 
have a continuous symmetry of the order parameter but 
lack "particle- hole" symmetryH), all or part of the disor- 
dered phase shows features which are usually characteris- 
tic of a critical point. More precisely, correlations in time 
decay algebraically and thus various susceptibilities may 
actually diverge. This behavior is due to Griffiths-McCoy 
singularitiegS^El, which arise from rare clusters which are 
more strongly coupled than the average, and it has be- 
come clear that their effect is much more pronounced 
near a quantum transition than near a classical critical 
point. 

One of the simplest models exhibiting the characteris- 
tic features of a quantum phase transition is the random 
transverse-field Ising chain, defined by the Hamiltonian 



L-l 



H 



L 



(1) 



where the {cr"} are Pauli spin matrices at site i and the 
interactions Ji and the transverse fields hi are random 
variables with distributions 7r(J) and respectively. 

A lot of results on the critical and off-critical properties 
of this model have been obtained, both analytically and 



numerically. The ground state properties of the Hamil- 
tonian in Eq. (|^) are closely related to a two-dimensional 
classical Ising model where the disorder is perfectly cor- 
related along one direction, the latter model first being 
studied by McCoy and WuEI. Subsequentljf. the quantum 
model was studied by Shankar and Murthyrl, and recently 
the critical properties have been worked out in great de- 
tail by D. S. FisJaer, using a real space renormalization 
group approachHE3. The quant una_piO|del, Eq. (|l|), has 
also been investigated numericallyt3~tj using a mapping 
to free fermions by means of a Jordan- Wigner transfor- 
mation. 

We now briefly summarize some of the surprising fea- 
tures of this model. Distributions of equal time corre- 
lation functions are found to be very broad which leads 
tea: (i) differjeat critical exponents for the divergence of 
the "typical" El and average correlation lengths, and (ii) 
the typical equal time correlation function at criticality 
falls off as a stretched exponential function of distance, 
quite different from the power law variation of the aver- 
age. Recently, it has been shown that these main results 
are not restricted to one dimension, but also seem t.n..hfpld 
in the two-dimensional random Ising ferromagnettStZI. 

A number of results for dynamics have also been found. 
For example, atj-the critical point the dynamical expo- 
nent z is infinitcu. Away from the critical point the dis- 
tribution of local (imaswiary) time dependent correlation 
functions is very broadO. The average varies as a (contin- 
uously vatving) power of imaginary time r involving an 
exponenttS z'{S), where 5 is the deviation from critical- 
ity. By contrast, the typical correlation function varies as 
a stretched exponential function of time. Quite detailed 
information on the whole distribution of time dependent 
correlation |£unctions in the paramagnetic phase has also 
been foundll3. At the critical point the average correla- 
tion function is found to decay with an inverse power of 
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the log of the timeO, corresponding to the result z = oo 
mentioned above. 

However, the distribution of time dependent correla- 
tion functions at the critical point has not yet been de- 
termined and is the focus of this study. We find that 
the distribution is very broad and can be expressed in 
terms of a single (logarithmic) scaling variable, defined 
in Eq. ^ below. This result implies that there is a con- 
tinuous range of exponents a, somewhat analogous to 
(but alsoj udth important differences from) "multifractal" 
behaviorEj'E3. 

We present our results in the following section, and we 
conclude with a discussion in Sec. [II. 



II. RESULTS 

Throughout the paper we assume the following rectan- 
gular distributions for the couplings Ji and the transverse 
fields hi 
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for < h < ho 
otherwise, 
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which are characterized by a single control parameter ho. 
The system possesses a critical point at S = [In J]av — 
[ln/i]av = 0, i.e. ho = 1, at which the distributions of the 
bonds and fields are equal. The lattice size is L and we 
impose free rather thaji .the more conventional periodic 
boundary conditionsMtJ. 

For the numerical work we make use of thej-mafiping of 
Hamiltonian (|l|) onto a model of free fermionsO^E^. Since 
the transformation has been used in previous work; .we. 
only give a brief summary here and refer to Refs.liililj'Lj 
for further details. For free boundary conditions, which 
we shall assume . hare , the most convenient representation 
is given in Refs.tliEj, necessitating only the diagonaliza- 
tion of a 2L x 2L real, tridiagonal matrix. The spin oper- 
ators occurring in the expectation value of the correlation 
functions can then be expressed as a product of fcrmion 
operators, which is evaluated using Wick's theorem. The 
resulting Pfafhan is then given by the square root of the 
determinant of a matrix, where the matrix elements can 
be calculated from the eigenvectors and eigenvalues of 
the Hamiltonian in the free fermion representation. The 
imaginary time correlation functions are always positive, 
so there is no ambiguity in sign when taking the square 
root. For L < 128 we average over 30000 realizations 
of the disorder, while for the largest size, L = 256, we 
average over 10000 realizations. 

We calculate the probability distribution P(lnS'(r)) of 
the single site imaginary time correlation function 



Suir) = (af(r)af (0)) 



(4) 



at the critical point, i.e. S — 0. For convenience, we will 
denote S'ii(r) by S{t) from now on. To obtain better 
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FIG. 1. The probability distribution, Ps, of -ln5'(r) for 
L — 128 and different values of the imaginary time r at the 
critical point 5 = 0. The data is averaged over 30000 samples. 



statistics, we determine the correlation function for ev- 
ery second site from L/4 to L/2, making a total of L/8 
sites. All sites are far from the boundary so we do not ex- 
pect the results to be affected by boundary effects. The 
average correlation function is then given by 



5'av(T) 



IS 



(5) 



where [• • -Jav denotes an average over samples. 

Since we expect strong finite size effects at the critical 
point, we calculated data for different system sizes. Data 
for the distribution of — In S{t) for L = 128 and L = 256 
is shown in Figs. |^ and||. One observes that the distribu- 
tions are broad and that for larger times the probability 
distribution gains more weight in the tail, indicating that 
correlations decrease for larger times, as expected. 

Since the curves for fixed L and different times r ap- 
pear to be shifted by a roughly constant amount on the 
(logarithmic) x-axis, we attempt a scaling plot of the data 
with the parameter free scaling variable 
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Note that in Ref. |l^, which investigated the dynamics in 
the paramagnetic phase, the scaling variable was found 
to be — lnS'(r)/T^/^, where ^ is expected to diverge at 
the critical point. Hence a in Eq. (|6|) is a natural scaling 
variable at the critical point. 

The corresponding scaling plots for P{a) against a are 
shown in Figs. || and ^ for L = 128 and L = 256. One 
sees that the data collapse is good for a not too large 
and that the range of a where scaling works increases 
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FIG. 2. The probability distribution of -ln5'(r) for 
L = 256 at the critical point 5 = 0. The disorder average 
is over 10000 samples. 
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FIG. 3. Scaling plot of the probability distribution in Fig. |l| 
(L = 128). The scaling variable a is that given in Eq. (^. For 
larger values of cx systematic deviations from scaling occur. 
This comes from data for small times r which is presumably 
not in the scaling region. 
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FIG. 4. Scaling plot of the probability distribution in Fig. 
^ (L = 256). The scaling variable a is that given in Eq. (|^. 
Note, by comparison with Fig. ^ that the range where the 
data scale well increases with increasing system size. 



with increasing L. The data for large a where scaling 
breaks down corresponds to small r, and it is reasonable 
to expect deviations from scaling in this region. 

Since the data is only a function of the scaling variable 
a, it follows that typically the correlation function falls 
off with a power of r. We shall now see that the precise 
value of the power depends in detail on how the typical 
correlation function is defined. For example, if we define 
"typical" to be the exponential of the average of the log, 



Se^vlogir) = exp([lnS'(T)]av) : 



we obtain 
[In^(r)], 



/ P{a)alnT da 
Jo 
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-(a)lnr, (8) 



which yields, the algebraic decay 
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where (. . .) denotes an average with respect to the dis- 
tribution P{a). From our data we get (a) ~ 0.7. 

On the other hand if we define "typical" to be the 
median of the distribution, then one easily sees that 
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median 
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where amcd is the median of P{a), i.e. it is defined im- 
plicitly by 



P{a) da. 



(11) 
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FIG. 5. An enlarged plot of the data in Fig. ^ (L=256). 
The full line is ~ a-°-^'^ . 



Any reasonable definition of "typical" will give a power 
law, in contrast to the average which has a much slower 
logarithmic variation, which we discuss next. 

Contributions to the average correlation function can 
come both frampthe scaling function and from non-scaling 
contributionalJH. The scaling part comes from the small 
a part of the scaling function where there is an upturn 
in the data, similar to that found for the scaUng func- 
tion of the distribution of the static spin-spin_correla- 
tions C{r) — {af{Q)af^^{0)) at the critical pointllfl. If we 
assume an algebraic relation P{a) ~ for small a, 
we can calculate the scaling contribution to the average 
correlation 5av(T) from 



P{a)S{T) da. 



(12) 



Noting that S{t) — exp(— alnr) one obtains 

/"OO 

5'av(''') ~ / a"'*' exp(— alnr) da 
Jq 



(Inr) 



-(l-A) 



(13) 



This agrees with the results of Rieger and IgloiEl who 
found 



of 0.67 is reasonably close to the value of A ~ 0.618 cal- 
culated above. Note again that there may be additional 
non-scaling contributions to the average correlation func- 
tion as in Rcf. nol 



III. CONCLUSIONS 

We have studied numerically the distribution of local 
(on-site) correlations in imaginary time for the random 
transverse- field Ising chain at the critical point. The dis- 
tribution was found to be logarithmically broad and the 
scaling variable a = log S{t) / log r was established. This 
means that while the correlations typically decay with a 
power of r, there is a range of exponents a with a distri- 
bution P{a). The small a part of the scaling function, 
which dominates the average correlations, can be fitted 
by P{a) ~ a""'^^, which gives (close to) the correct expo- 



nent in Eq. (13) for the logarithmic decay of the average 
correlation function. Note that all positive moments are 
determined by the small a region and so fall off with the 
samecB decay given in Eq. (|l3|) . 

The behavior of the distribution of S{t) found, here 
is somewhat analogous to "multifractal" behavioiO pre- 
dicted for the decay of spatial correlations in the jdnssical 
two-dimensional Potts model at the critical pointE], since 
both have a distribution of scaling exponents. Beyond 
that, however, there are significant differences. Whereas 
in our case the probability of having a scaling exponent 
a is P{a), which does not explicitly depend on r, for 
the corresponding multifractal behavior, the probability 
would be r"-'"'"^. As a result, for multifractal behavior, 
the typical correlation function has an exponent amin, 
the value of a at the minimum of f{oL), whereas here we 
find that the exponent for the typical correlation function 
depends on exactly how "typical" is defined. Averages of 
the n-th moment of the correlation function for positive 
n are also quite different. In our case, for all n > 0, the 
behavior is dominated by the small a region oi-the distri- 
bution and the exponent is independent of ?tE3, whereas 
for multifractal behavior, the m prnen ts depend on n in a 
non-trivial way and are given bylijnJ the Legendre trans- 
form of f{a). It would be interesting to see if there are 
other systems which have a "multiscaling" behavior of 
the type found here. 



IV. ACKNOWLEDGMENTS 



5av(r)^ (lnr)-2.„^ (^4) 

(where x„i — (1 — 0/2) ~ 0.191 is the bulk magnetization 
exponent with (p = [1 + \/5)/2), provided A = 1 — 2a;,„ ~ 
0.618. 

To check this we show in Fig. || an enlarged plot of the 
data in Fig. H (L = 256) for small a, together with the 
function ~ a , which is the best power law fit to the 
data. The curve fits the data fairly well and the exponent 
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